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NEW VERSION SUMMARY 



Title of new version: TOPAZO 2.0 
Catalogue number: 

Program obtainable from: CPC Program Library, Queen's University of 
Belfast, N. Ireland (see application form in this issue) 

Reference to original program: TOPAZO; Cat. no.: ACNT; Ref. in CPC: 76 
(1993) 328 

Authors of original program: Guido Montagna, Oreste Nicrosini, Giampiero 
Passarino, Fulvio Piccinini and Roberto Pittau 

The new version supersedes the original program 

Licensing provisions: none 

Computer for which the new version is designed: VAX, HP- APOLLO 7000; 
Installation: INFN, Sezione di Pavia, via A. Bassi 6, 27100 Pavia, and Sezione 
di Torino, via P. Giuria 1, 10125 Turin, Italy 

Operating system under which the new version has been tested: VMS, UNIX 

Programming language used in the new version: FORTRAN 77 

Memory required to execute with typical data 

300-600 kbyte as evaluator of observables in seven energy points 

No. of bits in a word: 32 

No. of processors used: 1 

The code has not been vectorized 

Subprograms used: NAGLIB [1] 

No. of lines in distributed program, including test data, etc.: 12343 
Correspondence to: 

GIAMPIERO@TO.INFN.IT, NICROSINI@PV.INFN.IT 
VAXTO : : GIAMPIERO , VAXP V: : NICROSINI 

Keywords: e + e~ annihilation, Bhabha scattering, LEP, Z° resonance, elec- 
troweak, extrapolated and realistic experimental set-up, QCD corrections, 
QED corrections, pure weak corrections, radiative corrections, Minimal Stan- 
dard Model, de-convoluted and realistic observables, 
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theoretical uncertainties. 
Nature of physical problem 

An accurate theoretical description of e + e~ annihilation processes and of 
Bhabha scattering at the Z° resonance is necessary in order to compare theo- 
retical cross sections and asymmetries with the experimental ones as measured 
by the LEP collaborations. In particular a realistic theoretical description, 
i.e. a description in which the effects of experimental cuts, such as maximum 
acollinearity, energy or invariant mass and angular acceptance of the outgo- 
ing fermions, are taken into account, allows the comparison of the Minimal 
Standard Model predictions with experimental raw data, i.e. data corrected 
for detector efficiency but not for acceptance. The program takes into account 
all the corrections, pure weak, QED and QCD, which allow for such a realistic 
theoretical description. 

Method of solution 

Same as in the original program. A detailed description of the theoretical 
formulation and of a sample of physical results obtained can be found in [2]. 

Reasons for the new version 

The new version gives the possibility of computing observables in an experi- 
mental set-up with different cuts on cross sections and asymmetries, according 
to the most recently published LEP data. Radiative corrections whose effect 
can become relevant in view of the present and foreseen experimental accuracy 
have been included. An option is added which allows an estimate to be made 
of the theoretical uncertainty associated with unknown higher-order radiative 
corrections of electroweak and QCD origin. 

Restrictions on the complexity of the problem 

The theoretical formulation is specifically worked out for energies around the 
Z° peak. Analytic formulas have been developed for an experimental set-up 
with symmetrical angular acceptance. Moreover the angular acceptance of 
the scattered antifermion has been assumed to be larger than the one of the 
scattered fermion. The prediction for Bhabha scattering is understood to be 
for the large-angle regime. 

Typical running time 

Dependent on the required experimental set-up. As evaluator of observables 
in seven energy points, between 10 (extrapolated set-up) and 270 (realistic 
set-up) CPU seconds for HP-APOLLO 7000, corresponding to 100-2500 CPU 
seconds for VAX 6410. As estimator of the theoretical uncertainty, it requires 
about 5400 CPU seconds for HP-APOLLO 7000, for a single energy point. 

Unusual features of the program 
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Subroutines from the library of mathematical subprograms NAGLIB for the 
numerical integrations are used in the program. 
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LONG WRITE-UP 



1 Introduction 

The high precision reached by the four LEP experiments has motivated sev- 
eral groups in assembling Monte Carlo programs or semi-analytical codes for 
giving accurate theoretical predictions relevant in the region around the Z° 
peak [1-6]. As a consequence we have by now several electroweak libraries for 
radiative corrections at LEP energies. This is a very important fact because, 
when precision physics is the main goal, continuous cross-checks are needed, 
especially at the moment when several new effects have been computed, due 
for instance to QCD corrections or large-m t behaviour 0. 

As it stands TOPAZO can be used to compute Z° parameters and de- 
convolved observables but also to obtain predictions for QED-dressed dis- 
tributions over a realistic set-up resembling the experimental raw data. The 
new version, TOPAZO 2.0, offers improvements both in the direction of comput- 
ing observables in a realistic set-up according to the selection criteria of the 
most recently published experimental data and in the direction of taking into 
account radiative corrections whose effect can become relevant in view of the 
experimental accuracy reached at present and foreseen for the near future. As 
far as the first kind of improvements is concerned, TOPAZO 2.0 allows one to 
put different cuts on the cross section and the forward-backward asymmetry 
in a given channel, thus matching the set-up adopted by some of the LEP 
collaborations. As far as the second kind of improvements is concerned, all 
recently computed electroweak and QCD corrections are implemented. More- 
over, further QED corrections whose size can become comparable with the 
foreseen experimental accuracy are included. 

For the time being, a particularly important item is to estimate the uncer- 
tainty intrinsic to the theoretical predictions. For instance as a consequence 
of the high experimental accuracy reached by the LEP collaborations, the the- 
oretical uncertainty affects the determination of the unknown parameters of 
the standard model |J. To this aim, a new feature has been implemented in 
TOPAZO, namely the possibility of running the evaluator over several theoretical 
options, related to different implementations of higher orders in the perturba- 
tive expansion. As a result, when choosing the proper branch in TOPAZO, the 
code returns an estimate of the theoretical uncertainty of electroweak and QCD 
origin. 

With respect to the previous version, the code has been rearranged in the 
form of SUBROUTINE, in order to render more viable the use of the program for 
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aims not planned by the TOPAZO package itself. 

2 The most important new features 

In the following we list and shortly comment the most important modifications 
of TOPAZO. They can be classified into three classes. 

TECHNICAL MODIFICATIONS 

• The present version allows the experimental cuts on cross sections and 
asymmetries in a given leptonic channel to be treated separately, accord- 
ing to the most recently adopted experimental strategy. To this aim the 
routines OBSERVABLES and OBSCUT have been updated. For instance the 
present version allows the computation of an extrapolated cross section 
and a cut asymmetry in the same run. 

• The new flag OTHERR has been included. For OTHERR = Y the evaluator 
runs over several (2 7 ) theoretical options, which reflect different (but not 
antithetic) implementations of radiative corrections, leading to a different 
treatment of missing higher-order terms. For OTHERR = Y the program 
returns a separate estimate of the weak and QCD uncertainty on the one 
hand and of the electromagnetic one on the other hand. 

• The code has been prepared in the following form of SUBROUTINE: 
SUBROUTINE TOPAZO (NRST , TRS , TZM , TTQM , THM , TALS) 

where 

NRST: the number of centre of mass energies 

TRS: the values of centre of mass energy (GeV) 

TZM: Z° mass (GeV) 

TTQM: top-quark mass (GeV) 

THM: Higgs boson mass (GeV) 

TALS: a s (M z ) 

are the basic input parameters. Moreover, the SUBROUTINE TINPUT is 
provided in order to allow the user to supply the experimental cuts for a 
given Z° decay channel and to select among different theoretical options. 
The meaning of the latter input parameters/flags is summarized at the 
beginning of SUBROUTINE TINPUT. In particular, as input/output facility, 
the flag OMON allows to write down the results for the de-convoluted 
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and realistic observables (OMON = Y), or to store them into appropriate 
common blocks /TTH , /TPO, /TPWTH (OMON = N). This second choice is 
supplied for fitting purposes. 



THEORETICAL IMPROVEMENTS: 
PURE WEAK AND QCD CORRECTIONS 

An important fact to realize is that there were several new calculations of 
radiative corrections in the last couple of years: their effect must be under- 
stood, classified descriptively, systematized and codified. Here it will only be 
possible to present a brief summary list of some of these effects: 

• By using a dispersion relation the value of a(M z ) obtained from experi- 
mental data is now a(M z ) = (128.87 ± 0.12)- 1 @. 



Two-loop heavy top effects in p,gv,A{b) for arbitrary values of m H (10 



0(aa s ) final state radiation ITT 



0(aa s ) corrections to vector boson self-energies [12 



Complete 0(a 2 s jjt) corrections to Ta, including NNLO terms for the 

z 

running b-quark mass |13 |. 



• Complete 0(a 2 ) corrections to Y z and 0(a^) corrections to Th |[3 

• 0(a s GFTn 2 ) corrections to Tf,, the FJRT effect [0 . 

• The 0(aa 2 ) correction to p |TjJ. 

To illustrate the numerical relevance of those effects included in our elec- 
troweak library we can compare the peak cross sections and the forward- 
backward asymmetries for p, r and hadronic channels in a fully extrapolated 
set-up and for the same input parameters as used in (see table 1). 
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Results 


TOPAZO 


TOPAZO 2.0 




1.4912 


1.4900 


O t 


1.4836 


1.4824 


Chad 


30.615 


30.591 




0.0047 


0.0039 


A fb (t) 


0.0047 


0.0039 



Table 1: The effect of the new weak and QCD corrections implemented in 
TOPAZO 2.0 for E c . m , = 91.222 GeV, M z = 91.175 GeV, m t = 200 GeV, 
m H = 250 GeV and a s = 0.125 



THEORETICAL IMPROVEMENTS: QED CORRECTIONS 

• The exact 0(a) initial-final state interference (including hard brems- 
strahlung) has been added for the leptonic channels as follows. The soft 
and virtual contribution has been included as a 1-dimensional spectrum 



(da/dcosfi, $ being the fermion scattering angle), taken from |16| and 
numerically integrated over the forward and backward hemispheres to 
build the correction to the cross section and the asymmetry. For the 
hard bremmstrahlung contribution, the 5-dimensional spectrum quoted 
in JT7J has been worked out analytically to obtain a twofold distribution 
[da jd cos ddE^, E 1 being the photon energy) and then numerically inte- 
grated as above. The analytical work performed on the hard photon part 
allows the computation of the interference correction with no substantial 
change in CPU time. The inclusion of the interference contribution is 
controlled, for each energy point, by a new flag 0NIF. For ONIF = Y, the 
new subroutine IFINT is called and the correction to the cross section 
and the asymmetry is returned after calling FUNSUBS and FUNSUBH for the 
soft and hard contribution, respectively. It should be pointed out that 
the treatment of the intial-final state interference is exact for s— channel 
annihilation but approximate for full Bhabha scattering. 

The 0(a 2 ) leptonic and hadronic contribution from initial-state radia- 



tion according to [18j has been included as follows. The soft and virtual 
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contributions are computed using the analytical formula given in p 



For the hard contribution, the 1-dimensional spectrum of |T8[ has been 
integrated numerically. The total relative contribution to the cross sec- 
tions in the hadronic and leptonic channels is returned for each energy 
point (for ONP = Y) after a call to subroutine PAIRS, which in turn calls 
the subroutines FUNSUBE, FUNSUBMU and FUNSUBHP, for e, fi and h 0(a 2 ) 
hard contributions respectively. 
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3 Test Run Output 

The typical calculations that can be performed with the new version of the 
program are illustrated in the following example. 



RESIDUAL WEAK CORRECTIONS ARE COMPUTED AT THE PEAK 
PAIR PRODUCTION EFFECT IS INCLUDED FOR EACH ENERGY 

I-F STATE INTERFERENCE IS INCLUDED FOR EACH ENERGY FOR LEPTONS ONLY 
S+T CHANNEL FOR ELECTRONS ARE COMPUTED 
NO CUT ON THE REDUCED ENERGY IS ASSUMED 

FULL PHASE SPACE FOR FINAL STATE EM RADIATION IS ASSUMED 
MINIMUM ANGLE OF MU(-) (DEG) : 0.00000E+00 
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MINIMUM ANGLE OF MU(+) (DEG) : 

MAX . ACOLLINEARITY ANGLE FOR MUONS (DEG) : 

MINIMUM ANGLE OF TAU(-) (DEG) : 

MINIMUM ANGLE OF TAU(+) (DEG) : 

MAX . ACOLLINEARITY ANGLE FOR TAUS (DEG) : 

ENERGY THRESHOLD FOR ELECTRONS (GEV) : 

MINIMUM ANGLE OF E(-) (DEG) : 

MINIMUM ANGLE OF E(+) (DEG) : 

MAX. ACOLLINEARITY ANGLE FOR ELECTRONS (DEG) 



O.OOOOOE+OO 
0.18000E+03 
O.OOOOOE+OO 
. OOOOOE+OO 
0.18000E+03 
O.IOOOOE+Ol 
0.40000E+02 
O.OOOOOE+OO 
0.25000E+02 



CURRENT VALUES FOR THE PARAMETERS ARE: 

Z MASS (GEV) = 0.91189E+02 TOP MASS 

HIGGS MASS (GEV) = 0.30000E+03 ALPHA.S 



(GEV) = 



. 17400E+03 
0. 12400E+00 



W MASS (GEV) 
NU 

ELECTRON 
MUON 
TAU 
UP 

DOWN (STRANGE) 

CHARM 

BOTTOM 



0.802951E+02 
0.167156E+00 
0.839136E-01 
0.839129E-01 
0.837222E-01 
0.300439E+00 
0.383139E+00 
. 300387E+00 
0.376022E+00 



SIN~2(E) 

SIN~2(B) 

A_FB(L) EFF. 

A_LR EFF. 

TOTAL WIDTH (GEV) 

G_H/G_E 

SIGMAO.H (NB) 

G(B)/G(HAD) 

A_FB(B) 

HADRONIC WIDTH (GEV) = 

A~POL_FB(B) 

A_FB(C) 

G(C)/G(HAD) 



0.232208E+00 
0.233490E+00 
0.150042E-01 
0.141608E+00 
0.249613E+01 
0.207727E+02 
0.414425E+02 
0.215718E+00 
0.991510E-01 
0.174312E+01 
0.934252E+00 
0.706791E-01 
0.172328E+00 
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DECONVOLUTED FB-ASYMMETRIES ARE: 
Z-EXCHANGE ONLY 



ELECTRON = . 1500028E-01 

MUON = 0.1500028E-01 

TAU = 0.1497573E-01 

CHARM = 0.7063442E-01 

BOTTOM = 0.9875908E-01 



DECONVOLUTED FB-ASYMMETRIES ARE: 
COMPLETE 



ELECTRON 

MUON 

TAU 

CHARM 

BOTTOM 



0.1628186E-01 
0.1628186E-01 
0.1629466E-01 
0.6814782E-01 
0.9661819E-01 



OBSERVABLES ARE 



E_CM (GEV) 



0.91300E+02 



SIGMA (E) (NB) 
SIGMA (MU) (NB) 
SIGMA (TAU) (NB) 
SIGMA (HAD) (NB) 
A_FB(E) 
A_FB(MU) 
A_FB(TAU) 
A_FB(C) 
A_FB(B) 



0. 1195475E+01 +/- 
0. 1486515E+01 +/- 
0. 1479019E+01 +/- 
0.3051606E+02 
0.1421349E+00 +/- 
0.7492082E-02 +/- 
0.7519941E-02 +/- 
0.6317556E-01 
0.9483699E-01 



0.1745525E-03 
0.8608965E-04 
0.8288089E-04 

0.1666541E-03 
0.7530680E-05 
0.7543220E-05 
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